Shell evolution and nuclear forces 
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We present a quantitative study of the role played by different components characterizing the 
nucleon-nucleon interaction in the evolution of the nuclear shell structure. It is based on the spin- 
tensor decomposition of an effective two-body shell-model interaction and the subsequent study of 
effective single-particle energy variations in a series of isotopes or isotones. The technique allows to 
separate unambiguously contributions of the central, vector and tensor components of the realistic 
effective interaction. We show that while the global variation of the single-particle energies is due to 
the central component of the effective interaction, the characteristic behavior of spin-orbit partners, 
noticed recently, is mainly due to its tensor part. Based on the analysis of a well-fitted realistic 
interaction in sdpf shell model space, we analyze in detail the role played by the different terms in 
the formation and/or disappearance of N=16, N = 20 and N = 2S shell gaps in neutron-rich nuclei. 

PACS numbers: 21.10.Pc,21.10.Jx,21.60.Cs 
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The shell structure is a common feature of finite quan- 
tum systems. Amongst them, atomic nuclei represent 
unique objects characterized by the appearance of a spe- 
cific shell structure. In particular, the magic numbers 
which correspond to the shell closures, will change de- 
pending on the N/Z ratio, i.e. when we move from nuclei 
in the vicinity of the ^-stability line towards the particle 
driplines. This has attracted a lot of attention nowadays 
because an increasing number of nuclei far from stabil- 
ity have become accessible experimentally (e.g., [l[ and 
references therein). The hope to reach even more exotic 
nuclei demands for an improved modelization, i.e. in 
the context of nuclear astrophysics. Since the underlying 
shell structure determines nuclear properties in a major 
way, changes of nuclear shell closures and the mechanisms 
responsible for that should be much better understood. 

Recently, the role of different components of the 
nucleon-nucleon (NN) interaction in the evolution of the 
shell structure has been actively discussed. Based on the 
analysis of the origin of a shell closure at N = 16, Ot- 
suka et al Q have suggested that a central spin-isospin- 
exchange term, f(r)(a ■ <r)(r • r) of the NN interaction 
plays a decisive role in the shell formation. 

However, from a systematic analysis of heavier nuclei, 
another conjecture has been put forward, namely, the 
dominant role played by the tensor force Q. The evi- 
dence is based on the comparison of the position of ex- 
perimental one-particle or one-hole states in nuclei ad- 
jacent to semi-magic configurations with the so-called 
effective single-particle energies (ESPE's). Within the 
shell-model framework, the latter ESPE's are defined [4j 
as one-nucleon separation energies for an occupied or- 
bital (or extra binding gained by addition of a nucleon 
to an unoccupied orbital) evaluated from a Hamilto- 
nian containing nucleon single-particle energies (the bare 
single-particle energies with respect to a closed-shell core) 



plus the monopole part of the two-body residual interac- 
tion [!, i| , i.e. 
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where j denotes a set of single-particle quantum numbers 
(nlj) and p refers to a proton (w) or to a neutron (y), n p - 

are particle-number operators. Vff, are centroids of the 
two-body interaction defined as 
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where the total angular momentum of a two-body state 
J runs over all possible values. 

The monopole Hamiltonian represents a spherical 
mean field extracted from the interacting shell model. 
Its spherical single-particle states, or ESPE's, provide an 
important ingredient for the formation of shells and in- 
terplay between spherical configurations and deformation 
in nuclei. Large shell gaps obtained from a monopole 
Hamiltonian are a prerequisite to obtain certain magic 
numbers. A reduction of the spherical shell gaps may 
lead to formation of a deformed ground state, if the cor- 
relation energy of a given excited configuration and a 
decrease in the monopole part are large enough to make 
such an intruder excitation energetically favorable. 

For example, the ESPE of the fO/7/2 orbital at Z = 8, 
N = 20 is the difference between total energy obtained, 
using Eq. ([T]), for 28 in its ground state and 29 with 
an extra neutron in the O/7 /% state assuming normal fill- 
ing of the orbitals (normal filling is used throughout this 
work). Considering a series of isotopes or isotones, it is 
clear that ESPE's will experience a shift provided by the 



2 



monopole part of the proton- neutron matrix elements, 
mainly. The bigger the overlap of the proton and neu- 
tron radial wave functions and the higher the j-values 
of the orbitals considered will lead, in general, to more 
drastic changes. In the present study we take into ac- 
count the mass dependence of the two-body matrix ele- 
ments of the effective interaction according to the rule: 
V(A) = (A core /Af^V{A core ). 

From the analysis of the experimental data and the 
ESPE's it has been noticed 0, H| that systematically 

l%l>l%l' I%l>l%l. ( 3 ) 

where j > =1 + 1/2 and j < =1 — 1/2 are proton orbitals 
and j> = V + 1/2 and j< =1' —1/2 are neutron orbitals. 
Thus, an extra attraction is manifested between general- 
ized spin-orbit partners (proton j = 1 + 1/2 and neutron 
j' = I' — 1/2 with I ^ I' or vice versa). 

This remarkable property is in line with the analytic 
relation valid for a pure tensor force |3|], i.e. using the 
above notation, (2j> + 1)V^., + (2j< + 1)V7^., = 0. 
To strengthen this idea, Otsuka et al. [3j have compared 
changes of the ESPE's in Ca, Ni and Sb isotopes, as due 
to the tensor force only and estimating its strength as 
resulting from a (tt + /^-exchange potential with a cut-off 
at 0.7 fm, with available experimental data. 

This work has stimulated a large number of investiga- 
tions using mean- field approaches [§-[23| . It is worth not- 
ing that phenomenological interactions, such as Skyrme 
and Gogny force, most frequently used in mean-field cal- 
culations, did not include a tensor term [24]. Provided 
its importance, a tensor term should be introduced and 
the parameters re-adjusted, what up to now, is not satis- 
factorily reached yet (see, e.g. Ref. dHHBj]). 

However, importance of the tensor force within the 
shell model 0, 11[ is mainly demonstrated in an empir- 
ical way. It is evident, that the choice of the partic- 
ular cut-off that was used to fix the strength of the 
tensor force component plays a crucial role in obtain- 
ing quantitative result for shifts in the ESPE's as pre- 
sented in Fig. 4 of Ref. 3j. It is also well known that 
the NN interaction is subjected to a strong renormal- 
ization before it can be handled as an effective interac- 
tion in many-body calculations within a restricted model 
space 26]. It is not straightforward to trace how the 
tensor component will become renormalizcd amongst the 
other terms contributing to the NN interaction. More- 
over, many shell-model interactions having high descrip- 
tive and predictive power were obtained by a x 2_n t °f 
two-body matrix elements to reproduce known experi- 
mental levels for a wide ran ge o f nuclei studied within a 
given model space (e.g. [27J, |28j). Even the effective in- 
teractions, maximally preserving their microscopic origin 
(based on a G-matrix), need further phenomenological 
correction (see e.g., [g, [20]). There is strong indication 
that inclusion of three-nucleon forces can heal the mi- 
croscopically derived effective interaction, in particular, 
improve its monopole part (see Ref. [3(| and references 



therein for ab-initio studies). However, there are still no 
systematic calculations available up to date for many- 
nucleon systems either within the shell model, or within 
the density-functional approach. This is why the present 
study of the two-nucleon case is still of interest. 

In spite of the indirect evidence at a two-body level [|| , 
up to now, the role played by the tensor force is not 
well determined. For example, recent shell-model studies 
based on large-scale calculations using a realistic effective 
interaction in the heavy Sn nuclei region [3l| conclude on 
the absence of a characteristic effect expected to result 
from a tensor force component. 

In this Letter we present a quantitative study of the 
role played by different components of the effective inter- 
action. It is based on the spin-tensor decomposition of 
the two-body interaction, which involves tensors of rank 
0, 1 and 2 in spin and configuration space. The pro- 
cedure allows to separate the central, vector and tensor 
parts of the effective interaction. The monopole prop- 
erties of each component can be studied separately, elu- 
cidating unambiguously its role in the shell evolution. 
The method has already been applied in a similar con- 
text [H, [HI, however, the authors used different effec- 
tive interactions in smaller model spaces, concluding on 
a second-order tensor effect only. Contrary to these re- 
sults, we put into evidence an important first-order tensor 
effect in the present study. 

A spin-tensor decomposition of the two-particle inter- 
action has been known for many years [3J-|40( . In a given 
model space, a complete set of two-body matrix elements 
determines the properties of nuclei ranging within this 
space. For spin 1/2 fermions (nucleons), one can con- 
struct from their spin operators a complete set of linear 
operators in a two-particle spin space: 

S^ = l,s! 2 0) = [a 1 xa 2 f\si 1) =a 1 +a 2 , 
= [ ai x a 2 } {2) , SP = [ ai x a 2 f } , = o x - a 2 , 

By coupling the spin tensor operators with the corre- 
sponding rank tensors in the configuration space one can 
construct scalar interaction terms. The most general two- 
body interaction can then be written as 

V(1,2) = V = (S (k) -Q {k) )= £ (4) 

fe=0,l,2 fc=0,l,2 

Here, and represent the central and tensor 
parts of the effective NN interaction. The V^ k=1 ^ term 
contains the so-called symmetric (S^g) and antisymmet- 
ric (Sj^-5 6 ) spin-orbit operators 37], which we will de- 
note as LS and ALS, respectively. To obtain the ma- 
trix elements for the different multipole components in jj 
coupling, first, one transforms two-body matrix elements 
between normalized and antisymmetrized states from jj 
coupling to LS coupling in the standard way. The LS- 
coupled matrix elements of can be calculated from 
the LS coupled matrix elements of V as 
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where a = (n a ,l a ). Finally, starting from the LS cou- 
pled matrix elements of V^ k \ for each k, we arrive at a 
set of jj coupled matrix elements to be used for further 
investigation. It is important to note that for a given set 
of quantum numbers of two-body states, the matrix ele- 
ments of V are a sum of the matrix elements of its three 
components V^ k \ 

In addition, using projection operators, one can se- 
lect different components of the effective interaction that 
connect two-nucleon states with specific values of the to- 
tal spin S, isospin T and parity (— 1) L . Thus, based 
on the selection rules in L^-coupling, we can separate 
triplet-even (TE), triplet-odd (TO), singlet-even (SE) 
and singlet-odd (SO) channels of the central part, as well 
as the even and odd channels of the symmetric spin-orbit 
and tensor part. 

Next, we propose to study the evolution of the ESPE's 
in series of isotopes and/or isotones induced by a given 
effective interaction and of its different multipole compo- 
nents. The decomposition described above is applicable 
only when the model space contains all spin-orbit part- 
ners within a given oscillator shell. This limits the region 
of applicability to the lighter nuclei. However, many in- 
teresting observations can still be extracted. 

In the present paper, we explore the effective interac- 
tion in IsOdlpOf shell-model space that reproduces very 
well the properties of stable as well as nuclei further away 
from stability [29j. We start with the "classical" case of 
the N = 20 isotones, going from O to Ca. In Figure 
1(a), we show the evolution of the neutron ESPE's us- 
ing the realistic interaction [29], relative to the energy 
of the 0d 5 /2 orbital. When protons fill the Oefe^ orbital 
(from O to Si), the splitting between neutron O/7/2 and 
O/5/2 orbitals decreases. The opposite effect is observed 
when protons fill the Ocf/3/2 orbital (from S to Ca) i.e. 
the corresponding splitting increases. A similar but less 
pronounced behavior is noticed for the neutron and 
lPi/2 orbitals. In Figures l(b-d), we show the variation of 
the same ESPE's, this time caused by the k = 0, 1 and 
2 multipole components of the full interaction. We re- 
mark that the summed shift of the energy for the various 
orbitals, produced by each two-body component, equals 
the total shift as produced by the two-body part of a 
full effective interaction, thus demonstrating the additiv- 
ity. We included in each figure the bare single-particle 
energies for visibility. 

From Figure 1 one concludes that the central part of 
the effective interaction barely produces any relative dis- 
placement of spin-orbit partners (Fig. 1(b)). A not very 
pronounced and often opposite effect is induced by the 
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FIG. 1: Variation of the neutron ESPE's in N = 20 isotones 
from O to Ca calculated using the realistic interaction [29l ] 
and its components separately. 



k = 1 part (Fig. 1(c)). It is indeed the tensor part which 
is responsible for the variation of spin-orbit partners, in 
line with the observation made in Ref. Q. In Table I 
we summarize contributions of the central, vector and 
tensor terms in the spatial even and odd channels, sep- 
arately. This elucidates the origin of the evolution of a 
given single-particle energy splitting. The decrease of the 
splitting between 5(^0^3/2) and e^ls^n) by 2.57 MeV, 
going from 28 to 34 Si (Table I, column 2), turns out to 
result from the combined effect of the central part (1.87 
MeV), in particular, in its triplet-even channel, and the 
tensor part of the nuclear interaction (1.06 MeV). 

Similarly, the increase of the gap between the neutron 
0^3/2 and O/7/2 orbitals when going from 28 to 34 Si and 
onwards from 36 S to 40 Ca (columns 3 and 4 of Table I) 
is a joint effect of the central and tensor component of 
the effective interaction. This is an important manifes- 
tation of the tensor force in this region. Due to the fact 
that at N = 20 the above two neutron orbitals have (i) 
the same radial quantum number, and, (ii) a different 
spin-to-orbital orientation, a large and negative tensor 
contribution of —1.93 MeV results for the variation of 
the gap between the 0d 3 / 2 and O/7 / 2 orbitals when filling 
the OCZ3/2 orbital with protons (from 36 S to 40 Ca). This 
large tensor shift, however, is almost fully cancelled by 



4 



I 

TABLE I: Contribution of different spin-tensor operators to the energy splitting variations A(j,j') = ej —Eji in different regions: 
N = 20 (columns 2-4), Ca-isotopes (columns 5-6), N — 28 isotones (columns 7-8). 



Energy (u0d 3/2 , v\ s 1/2 ) (v0f 7/2 , ^0d 3/2 ) (^0/ 7/2 , ^0d 3/2 ) (7r0d 3/2 , 7r0d 5/2 ) (7r0d 3/2 , tt1s 1/2 ) (^lp 3/2 , ^0/ 7/2 ) (^lp 3 / 2 , v0f 7/2 ) 



gap 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


Filling 
orbital 


7r0d 5 / 2 
28 0^ 4 Si 


7r0d 5 / 2 
28 0^ 4 Si 


7r0d 3/2 


^0/ 7/2 
40 Ca ^ 48 Ca 


v0f 7 , 2 
40 Ca ^ 48 Ca 


7r0d 5/2 
36 0^ 42 Si 


7r0d 3/2 
44 S^ 48 Ca 


Total 


-2.57 


3.68 


0.21 


-2.33 


-3.156 


1.60 


1.81 


Central 


-1.87 


2.17 


1.99 


-0.21 


-1.58 


2.03 


1.31 


TE 


-1.58 


2.23 


2.48 


0.62 


-1.19 


2.03 


1.02 


TO 


-0.68 


-0.31 


-0.11 


-0.03 


0.25 


-0.25 


-0.14 


SE 


0.71 


-0.45 


0.01 


-0.50 


-0.57 


-0.02 


0.18 


SO 


-0.32 


0.70 


-0.39 


-0.30 


-0.07 


0.28 


0.25 


Vector 


0.36 


-0.45 


0.15 


0.61 


0.06 


0.23 


-0.18 


LS 


-0.05 


-0.10 


-0.16 


0.09 


-0.15 


0.11 


0.15 


even 


-0.12 


-0.06 


0.25 


0.60 


0.25 


0.22 


-0.27 


odd 


0.07 


-0.04 


-0.41 


-0.51 


-0.40 


-0.11 


0.41 


ALS 


0.41 


-0.35 


0.31 


0.52 


0.21 


0.12 


-0.33 


Tensor 


-1.06 


1.96 


-1.93 


-2.73 


-1.64 


-0.67 


0.68 


even 


-0.78 


1.31 


-1.28 


-1.59 


-0.96 


-0.43 


0.43 


odd 


-0.28 


0.66 


-0.65 


-1.14 


-0.68 


-0.24 


0.26 



the central contribution of 1.99 MeV. The combined ef- 
fect results in only a slight overall decrease of the N = 20 
shell gap from 40 Ca to 36 S and 34 Si, thereby preserving 
the semi-magic nature of the latter nuclei. At the same 
time, while filling the 0d 5 / 2 orbital with protons (from 
28 to Si), due to the change in the spin-to-orbital ori- 
entation with respect the proton 0g? 3 / 2 orbital, the ten- 
sor contribution remains large but changes its sign (1.96 
MeV). This enforces the central contribution (2.17 MeV) 
and results in a rapid decrease of the N = 20 shell gap 
below 34 Si which is at the origin of the so-called 'island 
of inversion' around 32 Mg (deformed ground state). 

The position of the 0d 3 / 2 orbital and the possible 
shell gaps between this orbital and either the Is 1/2, or 
O/7/2 orbital, plays an important role in the formation of 
N = 16 as a magic number at Z = 8. This was formerly 
ascribed to result from the spin-isospin exchange part of 
the central force component [2[ and sometime later to be 
due to mainly a pure tensor force [4l|. The present re- 
sults support the important role of both a central part 
(in its spin-isospin exchange channel) and a tensor part 
in changing the shell structure between O and Si. 

In Figure 2, we show the two-body contribution to the 
binding energy from the monopole part of the realistic 
interaction and its different components. For the analysis 
we choose the N = 19 isotones 27 , 33 Si, 35 S and 39 Ca 
with a neutron hole (relative to N = 20) in the Od^/2 state. 
As is seen, the global shift is due to the central part of 
the effective interaction when adding up to 12 protons to 
the 16 core. This contrasts with the results, presented 
before, in which we studied the local relative variations 
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FIG. 2: Two-body contribution to the binding energy of O, 
33 Si, 35 S and 39 Ca in the lowest state (one hole in the 0d 3 / 2 
single-particle orbital), relative to the binding energy of 16 O 
(N = 20 isotones) using the realistic interaction [291 ] and its 
components, separately. 

in the single-particle energy in which the tensor force 
component plays a major role. 

More evidence results from the single-proton holes in 
Ca-isotopes i.e. studying the K isotopes. There is a cross- 
ing of the lsi/ 2 and the 0d 3 / 2 orbitals when approaching 
48 Ca and, in addition, a lowering of the energy gap be- 
tween the 0d 3 /2 and the 0rj 5 / 2 orbitals when going from 
40 Ca to 48 Ca, confirmed experimentally 0, 51 EI- In 
Figure 3, we show the variation of the proton ESPE's in 
Ca-isotopes obtained from the same effective interaction, 
while in Table I we present a detailed analysis of the role 
of different components in the evolution of the gaps. It 
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FIG. 3: Variation of proton single-hole states in Ca-isotopes 
using the realistic interaction [291 ] . 
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FIG. 4: Variation of neutron single-particle states in N = 28 
isotones using the realistic interaction [29| |. 



is seen (columns 5) that the lowering of the gap between 
proton 0d 5 /2 and 0d 3 / 2 orbitals as neutrons fill the Off/2 
orbital is mainly due to the tensor force. However, it is 
the central part, combined with the contribution from 
the tensor force, which reduces the gap between proton 
0d 3 / 2 and lsj/2 orbitals when approaching 48 Ca (column 
6). 

Finally, we explore evolution of neutron ESPE's in 
N = 28 isotones, from O to Ca, as a function of pro- 
ton number. The characteristic trends of spin-orbit part- 
ners, generic for a pure tensor force, are well manifested 
when using the same realistic interaction (see Fig. 4): 
approaching neutron Q/7/2 and Of 5/2 spin-orbit partners, 
and likewise for the IP3/2 and lpi/2 spin-orbit partners, 
when filling the proton 0d 5 / 2 orbital. An opposite effect 
results when filling the 0d 3 /2 orbital and fingerprints the 
contribution from a tensor term. 

In Table I (columns 7 and 8), we analyze in detail the 
reduction of the N = 28 shell gap, i.e. the change in the 
neutron IP3/2-O/7/2 energy difference from 48 Ca to the 
lighter isotones. To start with, these two orbitals have 
different radial quantum numbers. Therefore, the radial 
overlap contributing to the ^oJjg 2 n/ 7 2 centroid is larger 
than the radial overlap contributing to the Vq£ lpa a 
centroid. As can be seen from the table, the contributions 



from the central and tensor terms are dominating. Since 
both the O/7/2 and IP3/2 orbital are 'spin- up' oriented 
(j> = / + 1/2), the tensor term contributes in a similar 
way to the energy shift when protons fill the 0g?3/2 orbital 

(K^ _2 nf' y and V^~ 2 \Z U are both positive). The same 
happens when protons fill the 0rJ 5 / 2 orbital (Vq*~^o/t 2 
and Vod~2ip" 2 are both negative) . The overall difference 
in sign is due to the different relative spin to orbital ori- 
entation of the neutron orbitals (both (j > v = l v + 1/2) 
relative to the proton orbitals (j 71 = V ± 1/2)). Due 
to the difference in absolute value of the centroids, in 
particular, due to different radial overlaps for v0f-n0d 
versus vlp -?r0d, the positive and negative tensor contri- 
butions to the energy centroid do not cancel. They result 
in a shift of 0.68 MeV and -0.67 MeV filling the 0d 3/2 
or 0g?5/2 orbital, respectively. Adding this tensor energy 
shift to the central plus vector energy shift results in a 
reduction of the N = 28 shell gap going from 48 Ca to 44 S 
and from 42 Si to 36 O. This situation contrasts the N — 20 
shell gap evolution discussed before. 

These examples illustrate that in the discussion of shell 
gap evolution, it is mandatory to take into account what 
particular orbitals are considered. Both the central and 
tensor term represent important ingredients, together 
with the magnitude of the radial overlaps involved. 

To summarize, we have proposed a quantitative study 
of the shell structure evolution in series of isotopes or iso- 
tones, based on a spin-tensor decomposition of the two- 
body matrix elements. The method has allowed us to 
clarify the role played by the different terms of the ef- 
fective interaction in the variation of the single-particle 
energy of different orbitals. 

Based on the analysis of the best realistic interaction 
in the IsOdlpOf shell model space (29j, we show that 
the evolution of the N = 16, TV = 20 and N = 28 shell 
gaps is a combined effect of different spin-tensor terms, 
of which the central term in its triplet-even channel and 
the tensor term are of overwhelming importance. This 
conclusion partially supports the results of Ref. [12, HH 
regarding the importance of the triplet even channel but 
evidences the crucial role of the first-order tensor term 
as conjectured in Refs. [!, Q. The tensor term plays a 
dominant role, with increasing role of the vector term in 
the single-particle energy difference for spin-orbit part- 
ners. However, from the examples discussed here, one 
cannot assign unambiguously a dominating role to the 
tensor mechanism in cases when no explicit spin-orbit 
partners are considered. For example, the increase in en- 
ergy splitting between proton 0hn/ 2 and 0(77/2 when fill- 
ing neutron 0hn/ 2 in heavy Sb isotopes, or in the energy 
splitting between neutron Ohn/2 and 0(77/2 when filling 
proton 0<?9/2 in N = 51 isotones, discussed in Ref. [H, 
may be a result of different parts of the effective interac- 
tion. To clarify the observed situation, a corresponding 
quantitative analysis should be performed. 

The decomposition is a suitable tool only for the model 
spaces when all spin-orbit partners are present. This does 
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not allow, at the present moment, to analyze heavier nu- 
clei and check the hypothesis of the tensor force action 
in heavy systems, until a realistic determination of the 
position of relevant spin-orbit partners is established. As 
can be seen from the present work, the vector term of 
the effective interaction typically counter-acts the tensor 
term. The increasing role of the vector term could form 
a plausible scenario for a reduction of the tensor effect in 
heavy nuclei and explain the results obtained from shell- 
model studies of Ref. [3l[ . This requires more data on key 
heavy nuclei and the availability of extremely large-scale 
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